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I ORLWORD 

\n  experimental.  automatic  seamount  siiaiT  technique  for  satellite  rad.u  .iltmi.li>  has  ;  ecu 
developed  .itn!  evaluated  during  the  past  two  years.  It  was  based  on  an  cmpiiic.d  loriimla  in¬ 
ferring  the  seamount  slope  from  the  maximal  vertical  deflection  values  on  the  sea  suriu.c. 
Adjustment  of  the  maxima!  geoid  elevation  calculated  from  an  idealized,  theoretical  v.unount  model 
to  the  maximal  avoid  elevation  observed  on  the  altimetry  track  permitted  the  seamount  peal,  depth 
to  be  estimated.  Subsequently,  a  new  seamount  model  became  available.  Heine  much  more  de¬ 
tailed  than  the  one  just  mentioned,  it  made  possible  a  computer  study  invest  peat  me’  the  exact 
details  of  the  seamount  slope  ancle  estimation.  The  new  insights  into  the  nature  of  seamount 
parameter  estimation  then  led  to  the  development  of  a  second  detector  algorithm  documented  in 
the  present  report.  It  is  anticipated  that  the  new  algorithm  will  augment  rather  than  replace  the 
first. 

Dr.  B.  Zondek  of  the  Space  and  Surface  Systems  Division  contributed  the  matched  filter  as 
well  as  the  potential  theory  involved  in  the  interaction  of  the  seamount  and  its  root  with  the 
sea  surface.  Joint  Lilts  of  the  Physical  Sciences  Software  Branch  performed  the  computer  coding 
and  computer  program  checkout.  The  author  is  responsible  for  the  computer  program  formu¬ 
lation.  To  writ's  the  numerical  computer  program  checkout,  the  author  also  programmed  the 
algorithm  for  use  on  an  electronic  calculator  and  made  a  large  number  of  checkout  calculatioi  s. 

flic  work  documented  here  was  done  in  the  Space  and  Surface  Systems  Division  and  was 
funded  as  part  of  the  development  of  computer  programs  connected  with  the  evaluation  of  sea¬ 
mount  survey  techniques. 

Released  b\  : 

C).  I  .  HRAMOV  Head 

Strategic  Systems  Dep.iilmcnt 


NSWC  TR  H2-9I 


CON TUNIS 


IW 


INTRODUC  HON 


SPLCILIUAIION  Ol  1111  OM  -PARAML  H  R  SOI  11  ION 

1  OR  SLAMOUNT  SLOIM  .  3 

SAMI  LI  11  RADAR  ALUM  I  l'RV  DATA  FORMAL .  3 

MAK  111  l)  I  II  H  R .  5 

ROUGIINLSS  Din  (TOR .  <> 

SI  AMOUNT  LOCATOR .  ,s 

CALCULATION  OF  INITIAL  SLAMOUNT  PARAM1  T1  RS .  I.? 

OCl  AN  Dll’ I'll  ALGORITHM .  14 

I  SI  IMA  I  ION  Ol  PFAK.  DID  111 .  14 

ANALYSIS  Ol  DATA  DISPFRSION .  :<> 

AUTOMATIC  ANALYSIS  Ol  SATLLL1TL  ALTIMFTRY .  J‘> 


RFFLRLNUFS .  30 

APPFND1XFS 

A  TRIAL  DATA  LOR  ISOSTATIC  (  OMPLNSA TION .  31 

B  TRIAL  DATA  LOR  A  TYPICAL  CASF  OL  GFNI  RAL  (  OMPLNSA  LION .  43 

C  I  RIAL  DATA  I  OR  ABS1  NCF  OF  ROOT .  55 


DISTRIBUTION 


v 


NSWCTR  82-91 


INTRODUCTION 


An  experimental  computer  program  lor  seamount  model  parameter  estimation1 
became  available  last  year,  it  is  an  automatic  routine  designed  tor  application  to  SI  ASA  1  -A 
satellite  radar  altimetry  and  is  capable  of  performing  seamount  surveys  with  a  minimum  of 
human  intervention.  1  his  computer  program  is  based  on  an  empirical  formula  that  infers  the 
seamount  slope  iiom  the  maximal  vertical  deflection  values  on  the  sea  surface  above  the  sea¬ 
mount.  It  estimate-,  the  seamount  peak  depth  in  the  process  of  adjusting  the  maximal  geoi.i 
elevation  calculated  from  an  idealized,  theoretical  seamount  model  to  the  maximal  geoid  elevation 
observed  above  the  seamount  on  the  satellite  track.  Seamount  peak  geographical  coordinates  are  a 
byproduct.  Upon  successful  completion  of  the  numerical  checkout,  the  new  computer  program  wa> 
subjected  to  a  preliminary  test  and  evaluation  involving  several  Sl.ASAT-A  altimetry  tracks  for 

each  of  a  number  of  geophysically  diverse  ocean  areas.  That  study  was  severely  limited  by  the 
vast  gaps  in  the  coverage  of  the  globe  by  the  SUASAT-A  tracks.  Nevertheless,  the  results  of  this 
first  attempt  at  generating  quantitative  data  measuring  the  performance  of  the  experimental  sea¬ 
mount  detector’  indicated  that  fi'is  computer  program  can  he  e.xjrecteil  to  become  a  useful 
device  for  satellite  altimetry  data  analysis. 

1'he  principal  weakness  of  the  just-mentioned  experimental  seamount  detector  is  not  the 
simplicity  of  the  seamount  model  employed  (this  being  capable  of  calculating  the  geoid  deflection 
resulting  tront  the  presence  of  the  seamount  on  the  seamount's  symmetry  axis  only  )  out  rather 
the  empirical  formula  that  relates  the  slope  of  the  seamount  to  that  of  the  sea  surface  above 
This  formula  had  been  devised  by  scrutinizing  a  number  of  til  OS-3  radar  altimetry  data  tracks 

passing  over  seamounts  in  the  New  1  ingland  seamount  province.  The  study  had  been  confined 

to  that  part  of  the  world  because  a  rather  dense  net  of  altimetry  tracks  and  reliable  bathymetry 

were  noth  available  there.  As  part  of  the  work  preparatory  to  the  design  of  an  improved  de¬ 
tector  algorithm,  the  validity  of  the  slope  angle  formula  was  reinvestigated3  under  more  general 
conditions  than  those  prevailing  in  the  Western  North  Atlantic.  Although  an  equally  useful  combi¬ 
nation  ol  reliable  bathymetry  and  sufficiently  dense  ground  tracks  still  could  not  be  found 
elsewhere  in  the  oceans,  the  renewed  slope  angle  study  became  feasible  because  of  the  com¬ 
pletion  ol  a  new  seamount  model,  l'he  new  model  permitted  simulation  of  the  geoid  shape  above 
the  sc  imoi  nt  seamount  root  system  as  the  combined  effect  of  a  collection  of  gravitating  disks, 
the  cental  coordinates,  dimensions,  and  densities  of  which  could  he  individually  specified.  Most 
importantly .  this  new  computer  routine  was  aide  to  reproduce  geoid  height  and  vertical  deflection 
above  (lie  seamount  along  specified  straight  surface  tracks  intersecting  the  seamount  symmetry 
axis  or  offset  from  it  fo  obtain  a  detailed  relationship  between  seamount  slope  and  maximal 
vertical  dcllcction.  tour  drastically  differing  values  of  ocean  depth,  five  seamount  peak  sub 
mergence  depths,  and  various  typical  slope  angles  were  now  assumed,  l  or  each  of  altogethei  <>U 
different  model  seamounts.  each  consisting  of  20  disks  plus  a  corresponding  number  of  discs 
icpresenling  the  toots,  the  maximal  signature  slopes  were  calculated  and  compared  with  the 
assumed  slope  angles.  Additionally,  the  postulated  widths  of  the  seamounts  at  then  base  Were 
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compared  to  the  associated  signature  widths.  The  result  was  that  the  original  empirical  relation¬ 
ship  between  seamount  slope  and  signature  slope  proved  to  be  inaccurate  under  liie  variety  of 
conditions  considered.  On  the  other  hand,  it  was  realized  that  the  ratio  ot  signature  width  to 
seamount  width  at  the  base  is  fairly  insensitive  to  changes  in  the  seamount  parameters  and  ocean 
depth,  lhe  signature  width  being  known  from  the  altimetry  data  track,  it  now  became  possible 
to  calculate  the  seamount  width  at  the  base,  prior  to  the  actual  seamount  parameter  estimation, 
very  simply  and  ratiier  accurately. 

It  then  appeared  advisable  that  a  second  seamount  detector  should  be  developed,  starting 
with  tile  new  relationship  connecting  the  signature  width  and  the  seamount  base  width  and 
subsequently  proceeding  to  the  seamount  model  parameter  estimation  in  the  course  ot  which  the 
predicted  geoid  elevation  would  be  matched  to  the  observed  geoid  height  above  the  seamount.  In 
spite  of  the  usefulness  of  the  seamount  disk  mode!  for  the  purpose  of  the  new  slope  angle 
study,  no  particukh  advantage  could  be  realized  by  incorporating  it  into  the  new  estimator  algo¬ 
rithm.  The  original  seamount  model  was  thus  retained.  In  its  major  outlines,  the  new  seamount 
detector  is  similar  to  the  experimental  computer  algorithm  described  in  Reference  1.  Specifically . 
it  features  the  digital  filter  necessary  to  recognize  and  enhance  the  seamount  signatures  concealed 
among  the  radar  altimetry  data,  the  basic  model4  containing  the  seamount  related  potential 
theory,  and  a  mathematical  estimator  for  the  physical  seamount  parameters.  The  latter  estimator 
adjusts  the  characteristic  properties  of  the  radar  altimetry  signature  predicted  from  the  theoretical 
model  to  the  corresponding  quantities  among  the  empirical  data.  The  present  and  the  experi¬ 
mental  seamount  detectors  differ  mainly  in  this  latter  computer  program  segment.  While  the 
previous  algorithm  performs  at  this  point  a  one-parameter  solution  for  the  seamount  base  width, 
the  present  version  of  the  seamount  detector  first  calculates  the  seamount  base  width  from  the 
signature  width  and  then  iteratively  varies  the  remaining  seamount  dimensions  (in  particular,  slope 
angle  and  peak  submergence  depth)  until  the  computed  maximal  geoid  elevation  above  the  sea¬ 
mount  matches  the  altimetry  signature  height  within  a  specified  margin,  subsequently  attributing 
the  resulting  peak  depth  and  slope  angle  to  the  actual  seamount.  Again,  because  of  the  imper¬ 
fections  inherent  in  our  present  knowledge  of  seamount  compensation,  this  estimation  is  done 
twice  per  detection,  separately  for  the  cases  of  perfect  isostasy  and  absence  of  root,  with  an 
option  to  prescribe  arbitrary  compensation  if  desired.  The  actual  depth  value  is  once  more  as¬ 
sumed  to  be  bracketed  by  the  depth  values  resulting  from  the  isostatieally  compensated  and 
uncompensated  cases. 

It  should  be  noted  that  a  change  was  specified  for  the  experimental  seamount  detector 
after  it  was  documented,'  largely  as  a  result  of  the  work  leading  to  the  present  algorithm.  This 
change  concerns  the  computer  program  segment  "Roughness  Detector”  (pages  <s  and  l)  of  Ref¬ 
erence  1),  in  particular,  the  definition  of  the  ”ina.xiimim-SNv-wiiuknv."  This  segment  of  the 
algorithm  lias  been  deleted,  it  was  replaced  by  a  new  version  that  is  identical  to  the  "Roughness 
Detector”  specified  below  as  part  of  the  present  "One-Parameter  Solution  for  Seamount  Slope.” 
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SPECIFICATION  OF  THL  ONE-PARAMETER  SOLUTION 
FOR  SEAMOUNT  SLOPE 


SATELLITE  RADAR  ALTIMETRY  DATA  FORMAT 


Flic  satellite  altimetry  uata  to  tie  processed  try  the  seamount  detector  can  he  recalled  Iron, 
permanent  storage  in  the  format  indicated  in  Tables  1  and  2  (Ci.  It.  West,  unpublished  d  Pa. 
D>7S).  Extract  from  the  data  file  header  the  time  interval.  At.  in  seconds  and  the  subsalellite 
velocity.  \s.  in  km  sec.  From  each  segment  of  the  altimetry  data  record,  obtain  for  cacti  uata 
point  (time  instant.  tr  in  seconds)  the  filtered  geoid  height.  Nj  =  .\(t,l.  in  meters,  the  deflection 
of  the  vertical.  5  =  6(1!).  in  arc  seconds,  the  latitude.  <£,  -  ght,).  in  radians,  and  the  longitude 
last.  ,\  =  ,\(  t, ).  in  radians.  Convert  ^  and  A,  to  decimal  fractions  of  degrees.  Note  that  the 
altimetry  data  serving  as  input  to  this  seamount  detector  are  not  the  raw  altimetry  data  u.icks 
but  result  themselves  from  a  filtering  process  explained  in  References  22  and  22  of  NSVvt  IK 
S1-2U0.1  Here  and  in  the  following  we  shall  refer  to  N|  and  6,  as  "filtered”  when  wishing 
emphasize  their  origin,  try  Kalman  smoothing,  from  the  raw  satellite  radar  altimetry,  \tiih 
reference  to  the  present  seamount  detector,  and  6,  will  usually  be  considered  "untiltereu”  and 
named  so  because  they  serve  as  input  data  to  the  matched  filter. 
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5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 


TABLE  1. 


Type 

1 

I 

1 

1 

R 

R 

A 

1 

R 

R 

R 

R 

R 

K 

R 

R 

R 

K 

K 


SEASAT-A  FILTERED  GEOID  HEIGHT  DATA  FILL  (TYPE  lGD) 
HEADER  RECORD  (COMMON  HEADER) 


Approximate  Range 
of  Significance 


Description 


XXXX 

XXXXX 

XX 

XXX 

XXXXX.XXX 

xxxxxx 

YYDDD 

XXX 

XXXX 

XX.XX 

XXXX 

XX.XX 

X.XX 

±xxx 

XXXX 

X.XX 

xx.x 

±.xxx 


Number  of  points  in  file 
Rev  number 

Starting  year  of  segment 
Starting  day  of  segment 
Starting  seconds  of  segment  (see) 

Time  interval  (sec) 

't  ear  and  Julian  day  (right  adjusted) 

Altimeter  mode 
Autocorrelation  distance  (km) 

Standard  deviation  of  data  (m) 

Standard  deviation  of  geoid  heights  (m) 

Standard  deviation  of  vertical  deflections  (are  see) 
RMS  of  filtered-raw  differences  ( m! 

Maximum  vertical  deflection  ( arc  see) 

Average  velocity  (km/sec) 

Antenna  distance  from  satellite  center  of  gravity  (in' 
Radar  instrument  delay  distance  equivalent  (cm) 
Time  correction  (sec) 

Sparc 
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TABLL  2.  SLASAM  FILTERED  t.EOID  MIGHT  DATA  I  il  l:  (TYPE  1  t.l)> 
DATA  RECORD  (COMMON  DATA) 


Word 

T  y  pc 

Approximate  Range 
of  Significance 

Description 

1 

R 

X.XXXXX 

Latitude  ( rad  i 

•> 

R 

WWW 

L  ongiludc  ( rad  ) 

1 

K 

•  w.\.  w 

< icoid  height  (  m) 

4 

R 

WWW 

Vertical  deflection  (aic  seel 

5 

R 

:X\\.W 

Raw  gcoid  height  t  m  i 

0 

R 

XXX 

(icoid  height  confidence  hound  tin) 

T 

R 

XX  XX 

\’e rticaJ  deflection  confidence  hound  (arc  seel 

,s 

K 

X  XX 

Orbit  uncertainty  ingeoid  heights  (ml 

9 

R 

XXX 

Orbit  uncertainty  in  vertical  deflections  (arc  sec) 

U) 

1 

>1  -->14 

f  lag  word 

1 1 

R 

XX  X 

Significant  wave  height  (tin 

I  2 

R 

XX.  X 

o  -  S\V|  1  ( m  i 

15 

R 

XX 

Automatic  gam  control  ( d  13 > 

14 

R 

XX 

a  -  At (dlD 

15 

R 

±  XX 

Tilt  SWT  1  correction  ( m  l 

In 

K 

x.x.x.x 

lilt  1  rad) 

17 

R 

XX.  X 

Ionospheric  correction  (cm) 

18 

R 

XXX 

Atmospheric  pressure  (mb) 

19 

R 

XXX.  X 

Dry  tropospheric  correction  (cm) 

1  J 

o 

R 

XX.  X 

Wet  tropospheric  correction  (cm) 

21 

R 

±xx.xx 

Tide  ( m  l 

1  A 

R 

x.xx 

Barotropic  correction  (m) 

25 

R 

XX 

Wind  speed  ( kn  I 

24 

R 

X..X 

Wind  direction  (  rad ) 

25 

R 

tX.X.X 

Sea  state  tSWli)  correction  (m) 

28 

R 

XX.  X 

Alternate  wet  tropospheric  correction  (cm) 

27 

R 

XX.  X 

Rain  rate  (  mm  hr) 

28 

R 

XXX.  X 

Steric  correction  (cm) 

29 

R 

Spare 

50 

R 

Sparc 
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M  A  I  CHI  D  FILTER 

l  he  male1  \!  filler  is  an  optimized,  digital.  high-pass  filler  permanentlv  stored  on  tile  ha.i. 
eomputer  svstem.  lls  application  to  the  Kalman  filtered  SI  ASA  I -A  altimeliv  data  base  for  ttie 

purposes  of  the  seamount  deteetor  ls  controlled  in  the  SI  ASA  I- \  revolution  numivt  the  devn.e 
set  number  associated  with  the  particular  revolution  inimhcr.  and  n\  the  fiitei  constant  Inc 
latter  is  related  to  the  cutofl  wavelength  \  ot  me  tiller  anu  the  length  I  oi  the  sampling 
interv  al. 

,7  I 

p  tail  —  i  , 

\  and  1  are  specified  in  terms  of  km.  \  tv  picul  value  for  SI- ASA  I  -  A  data  is  I  5.57.  >m 
tv  pica!  value  pairs  lor  and  p  are 

,v  =  50  km  .  p  -  u.:i: 

\.  =  10U  km  .  p  =  0.1047 

X,  =  150  km  .  p  =  0.00%  . 

Adopt  the  following  ternimologv  N  are  the  "unfiltered  geoid  height"  values,  anu  «>  are 
the  "un filtered  vertical  deflection"  values.  Application  of  the  matched  filter  to  \  and  X  will 

/\  \  a  'a  1 

result  in  the  N  and  6  data  tracks.  \  are  the  “filtered  vteoid  heights"  i  m  mi  and  A  a  re  t.ie 

i  i  i  ' 

“filtered  vertical  deflections"  (in  arc  sec).  I  he  index  “i"  is  related  to  the  time  value  i  associated 

1 

with  the  data  values. 

Note  that,  because  of  the  nature  of  the  filter  algorithm,  no  filtered  data  values  will  be  as¬ 
sociated  with  the  first  few  and  the  last  few  unfiltered  data  points.  Thus,  the  first  “data  window" 
(see  below  for  definition)  on  a  data  track  shall  start  with  the  first  available  Altered  .  :oid  height 
and  vertical  deflection  and  the  last  data  window  shall  end  with  tne  last  available  pair  of  filtered 
data  values,  further,  the  description  of  the  control  card  deck  and  data  cards  on  pages  0  ami 
of  Reference  1  arc  now  obsolete  and  must  be  replaced  b\  the  cards  listed  in  la'Me  3. 

TABLE  3.  CONTROL  CARDS  I  OR  MATCHED  FILTER 

JOB  C  ARD  ( 1*4.  I  I  10) 

U'COCNT  C  ARD 
A  IT  AC  1 1.  J(  1  .  ID  \B/ 

BIT  if  N.  BIII’SA.  .1(1.  Rl  \  =\.\\\.SN=r\L-=^.()=OC  ]  .1  (  =CC.RII()E  =  YYY.\  e.l’l  ( )  I  NO 


(1  ~  |S  o 

\\\\  .  f  our-digit  revolution  number  (Examples  0747.  1754.  els.  i 

N  l '»*•?■=■  .  Device  set  number  associated  with  SI  ASA  I  file  i-oDOWWI 

YYY  .  f  ractional  part  of  the  input  parameter  “rho" 
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I  hesc  control  cards  will  pioducc  till  il.iLi  iik-  "I)  \  I  HI \  \ \  WKI U»'»  )  \  NS"  idcntiiicu  m. 

I  Ik'  ID  specified  on  ilio  ion  card.  When  mtciuang  to  inn  tli  vjinmini  UcUkioi.  inis  til.  iv 
attached  as  I  Al’l  I. “  "I  C  — C  C  "  m  Ilk-  lourth  control  ^aru  iiioib  that  the  printed  output 
generated  hv  the  fillet  algorithm  kill  he  leprodiiceU  on  nikioikiic. 


ROLt.HNLSS  1)1  I  I  t  TOR 
Input  Data 


N,  Data  track  1 1  iltcic.l  v  •  s  :  in 


Data  track  Milk  ;. ..  v .  :.  :k 

,  l-  * 

JlL  SC‘l 

At 

Data  step  vviut'i.  Al 

M/v 

V  s 

Suosalellite  vchkitv 

km  see 

L 

l)ata  window  widtn  isp.citic!  nv 

i i nc  r  i 

kill 

deln 

1  hrcslloki  par.uneki  ispculkki  h 

\  itM'II 

m 

1 

Data  window  widtn  ’actor  (user 

UV.  llied  1 

Default  \alue  1=1.'  dimensionless 


Algorithm 


Consider  the  filtered  data  tracks  N  and  6  .  Assien  the  value  zero  to  ,dl  \  that  are  negative. 

ii-  i 

Specify  data  track  "windows."  The  number  of  the  data  points  within  each  window  is 


n  =  1M 


+  1 


(501) 


Himinate  all  data  track  portions  that,  after  application  of  the  matched  filter,  have  a  number  of 
points  that  is  less  than  1.5  to  3.U  times  the  number  n  of  data  points  within  the  window.  Start 
the  first  window  at  the  first  data  point  not  omitted  by  the  filter,  t  and  calculate 


Do  the  same  for  the  next  window.  N  .  \k+,  . Nk  +  n.  obtaining 


1 5  o :  i 


(303) 


(i 
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ami  tor  each  hollowing  window,  obtaining 


SN 


k>  > 


i  .'>04  i 


until  the  data  track  \  is  exhausted. 

I 

from  now  on.  consider  only  those  data  windows  lor  which  SN  '  deln.  Note  that  th.  latt.-i 
requirement  divides  the  new  data  track  (consisting  ot  the  SN  %  allies  I  into  "islands."  lea  aid, 
data  island,  calculate  the  center  of  gravity. 


SN  =  SN  (t  i  <ao5. 

1  1 

It,  SN, 

1<'<U'  2;  SN, 

where  the  summation  extends  over  the  island  under  consideration.  Note  that  tr0(,  will  generally 
not  exactly  coincide  with  a  t,  data  value.  Once  tC0(;  has  been  found,  substitute  for  it  the  t 
value  directly  preceding  it  or.  if  more  convenient,  the  t,  value  nearest  to  tC  O(; :  tt  (Ki  -  t*(Mi 
from  all  data  windows  associated  with  the  data  island,  select  then  that  window  that  is  charactei- 
i/ed  by  t£0(j  (the  first  data  point  of  which  has  the  time  value  t*ot;  )• 

Lx  tend  finally  the  thus-chosen  data  windows  (one  for  each  data  island)  to  the  left  In  m 
data  points  and  to  the  right  by  m  data  points  (i.e..  by  data  points  taken  from  among  tin 
just-calculated  SN(  values).  The  number  m  is 


Output  Data 


tetri 


Lor  each  of  the  latter,  extended  windows,  store  and  print  out 
t  t  ..  t  . 

M  +  I  n  +  n  -  1 

A  A  A 

N  N  .  N  , 

u  M +  •  m  +  n  -  1 

A  A  A 

5  6  .  5 

M  M  +  1  V  +  n  -  » 

SNV 

N  N  .  N  , 

6  6  .  5  . 

m  n  + 1  u  +  n  -  J 
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SEAMOUNT  LOCATOR 
Subroutine  "LXfU" 


Known  are  three  data  points,  i  =  1.  2.  3.  on  a  data  track  \  =  vt\  I.  Assume  that  these 
points  bracket  a  nu\muun  or  niinimuni  ot  the  \  values.  I  he  abscissa  vaiue.  \  .  associated  will) 

- 1  i' 

the  extremum  ot  \  is  then 
*  1 

M 

\  =  \.  -  — -  (4G1 l 

i'  1  2cd 


>:  -  >t 

t\:  - 

•V2 

bd  = 

(402  > 

y.r  -  >  I 

(X.,  - 

Xl 1* 

-  X, 

V:  - 

y , 

cd  = 

(403) 

'-a  - 

y_,  - 

>1 

Subroutine  “INTERPOL-I ” 

Known  are  three  data  points,  i 

=  1.  :.  3. 

on  a  data  track  v  =  vlx.l 
"1  '1 

.  hind,  for  a  specified 

abscissa  value.  x  .  the  ordinate  value 

p 

>p  =  ><V 

y  =  a  +  b  (x  - 

P  P 

V  + 

(404) 

a  =  y, 

(405) 

1  -  >1 

(x2  - 

x,)2  | 

b  =  - 
d 

y.)  -  y  i 

(x3  - 

Xl)2 

( 40(r ) 

,  :x2  - 

yr  - 

Yl 

1  i 

c  -  ~ 

.1 

(407) 

LI 

*3  -  *1 

>’j  - 

yi 

d  —  (X2  ■*  X  j  )  (X3 

-  *1 >  (Xj  - 

-  x2) 

(408) 
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Subroutine  ■  INTI  RPOL-2 


k  in 

uvn  are  tour  c 

ata  points,  i  =  1 

.  3.  4. 

on  a 

data  track  y(  =  y(x  ).  Calculate,  tor 

specified 

abscissa  value. 

x  .  the  ordinate 

value  y  = 

y(xe 

). 

V 

■  r 

-  c  +  It  x  -  > 

p 

i  >  +  g(xp  -  x,  )2 

+  h(.x  - 

p 

x,  )3 

(409i 

e  “ 

>  i 

(410) 

i 

>  :  -  y  1 

t.x:  -  x ,  )- 

(Xj  -  X 

f  =  — 

y.»  -  >  i 

y4  -  Vi 

IXj  -  x,  >2 
t  x4  -  X,  )2 

(x3  -  X 

(x4  -  X 

t3 

>3 

(41  1  i 

i 

\:  -  X, 

y  2  ~  y  i 

(Xj  -  X 

)3 

i 

A 

Xj  -  x, 

X4  -  X, 

y 3  -  yi 

>4  "  >1 

(Xj  -  X 

(x4  -  X 

I3 

)3 

(412) 

1 

X;  -  X, 

(Xj  -  X,)2 

y2  -  yi 

i 

h  =  - 

\ 

Xj  '  X, 

(Xj  -  X,)2 

yj  -  yt 

(413) 

\4  -  \  , 

<x4  -  X,  )2 

y4  -  yt 

x2  -  X, 

(Xj  "  X,  )2 

(x2  -  X 

,)3 

A= 

Xj  "  X, 

(Xj  -  X,  )2 

(Xj  -  X 

,)3 

(4 14  > 

x4  -  X, 

(X4  -  X,  )2 

(x4  -  X 

.)3 

Note  Concerning  the  Subroutines 

Optionally,  the  computer  programmer  may  implement  the  above  subroutines  by  any  of  the 
classical  interpolation  schemes,  if  desired,  the  polynomial  coefficients  can  be  evaluated  by  a  least- 
squares  fit  method.  Irrespective  of  the  method  chosen,  it  must  however  be  kept  in  mind  that 

rather  frequently  along  the  altimetry  track,  the  argument  values  (time  or  distance  -  depending  on 
the  interpretation  of  x.)  will  be  large  compared  to  the  time  required  for  a  seamount  overflight 

or  to  the  seamount  width,  implying  x  >  (x  -  x  ).  To  avoid  loss  of  significant  diaits  durine 

the  evaluation  of  the  determinants,  it  is  advisable  to  first  transform  the  x.  values  to  corresponding 

xj  values  so  that  the  origin  of  the  x'  system  coincides  with  Xj .  xj  =  (x(  -  x t  >.  implying  xj  =  0. 
1  he  particular  interpolation  method  chosen  may  then  be  executed.  Subsequently,  the  results  must 
be  transformed  back  into  the  x/y  coordinate  system.  Note  that  this  precaution  is  unnecessary  if 
the  subroutines  are  applied  in  the  form  specified  above. 
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Input  tor  Seamount  Locator  Algorithm 

1  or  each  data  window,  have  available 


Nil) 

i  i 

Cntiltered  geoid  height 

m 

6  1 1  ) 

1  1 

In  filtered  vertical  deflection 

arc 

t 

] 

1  ime  argument 

see 

X 

1 

L  ongitude  l  ast  from  altimetry  file 

dCt ! 

'r. 

Latitude  from  altimetry  tile 

d  Cl! 

Seamount  Locator  Algorithm 

Lor  each  of  the  specified  data  windows,  now  consider  the  unliltered  geoid  heights.  N  .  and 
the  unt'iltered  detlections  ol'  the  vertical.  6  Also  keep  in  mind  that  it  is  possible  to  associate 

with  the  discrete  number  sets  Mm  and  <\tM.  the  continuous  functions  Nltl  and  6<ti  representing 

the  geoid  height  and  vertical  deflection  in  the  “real  world"  or  simply  being  assumed  to  be 
suitable  functions  adapted  to  N  it  t  and  6^1  i  In  least-squares  fit  or  other  meaningful  methods, 
l  or  the  present  purpose,  it  is  sufficient  to  postulate  that  Ntt)  and  6m  may  be  realized  mathe¬ 
matically  if  needed.  Perform  the  following  calculations. 

1.  When  traversing  the  window  from  the  left  to  the  right  (proceeding  from  lower  values  ol 

t  to  hieher  ones),  expect  to  encounter  negative  6  values  and.  aimmi:  them,  an  absolute 
minimum  value.  Consider  the  data  point  5  (t  )  associated  with  this  minimum  value  plus 
the  data  point  to  the  left  and  the  data  point  to  the  right  of  it.  Apply  to  these  three 

data  points  the  subroutine  "1\1R"  to  find  the  abscissa  value  t^  belonging  to  the 

minimum  ol  the  curve  6(il  associated  with  the  just  specified  three  data  points. 

2.  In  the  same  manner,  find  the  tune  value’  t  associated  with  the  absolute  maximum 
value  of  6  It  )  that  will  occur  further  down  the  data  track  but  still  within  the  data 

I  i 

window. 

3.  Within  the  window,  the  data  track  \<t  )  mav  ne  expected  to  have  an  absolute  maximum. 
Lind  the  data  value  representing  this  maximum.  Also  find  the  neighboring  data  points 
to  the  left  and  to  the  right  ol  this  maximum.  Apply  subroutine  ‘l  X I  K  "  to  find  the 
time  value  tSM  of  the  maximum  of  Nltl  related  to  the  just-defined  three  data  points. 


NSYVCTR  82-91 


4.  Now  apple  subroutine  "IN  11  Rl'OL-l"  to  tSM  and  the  three  data  points  on  the  \ 
track  from  which  tSM  was  determined  in  step  3.  obtaining  \(tSN)i  -  N 

5.  Nest,  determine  tile  four  data  points  on  tire  Ni  track  that  are  nearest  to  (lie  time 
value  t  Apple  to  them  subroutine  'IN  T1  Rl'OL-2"  to  evaluate  the  \  value  associated 
with  t4  .  Nit.)  =  N . 

0.  further,  proceed  as  in  step  5  to  find  Ni  t  B )  =  \  . 

7.  Proceed  as  in  step  5  to  find  8(tA>  =  &A . 

8.  Proceed  as  in  step  6  to  find  6(  tu  )  =  8^. 

9.  Consider  now  the  data  tracks  N.  and  8.  in  the  neighborhood  of  lSM.  in  particular,  the 

time  value  t  =  t~M  just  to  the  left  of  (prior  to)  tSM  and  the  lime  value  t.  =  t’  to 
the  right  of  (following)  tSM.  Note  that  t~M  <  tSM  and  t*V)  >  t^.  In  the  radar  tiles. 
t~M  and  t*M  are  associated  with  the  values  of  longitude  Last  and  latitude  on  the 
subsatellite  track. 


'sM  ''‘SM 

and 

*SM 

f  +  .  \  + 

lS  M  AS  M 

and 

^SM 

Evaluate  the  longitude  Last  and  the  latitude  of  the  "signature  peak. 


^SM  ^SM  4 


M  ~  M 


fS  M  *SM 


ltSM  'sM  1 


(415i 


10. 


^SM  '  ^SM 

^SM  ^SM  +  ~ 7  (  tSM  “  lS  M  ' 

t  - 
lSM  lSM 


<4 1  o  i 


Note  that  above  the  data  were  postulated  to  conform  to  the  typical  seamount  altimetry 
pattern  illustrated  in  Figure  1.  If  the  data  fail,  altogether  or  in  any  essential  detail 
to  match  this  pattern,  reject  the  data  window. 


Seamount  Locator  Output 


For  each  data  window  successfully  processed,  store  and  list  the  following. 


1 1 
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P 


A 


URL  I.  TYPICAL  CiLOII)  1 1  IK.IIl  PATTI  RN  AND  VI  K  IK  AL  Dll  LI  Cl  ION 
ABOVL  A  SLAMOUNT 


1: 


J 
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N. 


!sm  \  m 


\s  M 


SM 


Seamount  longitude  and  latitude  are  to  he  specified  in  decimal  degrees  as  ncll  as  in  degree  minute 
-  second  format. 


CALCULATION  Ol  INITIAL  SEAMOUNT  PARAMETERS 


Input  Data 


N 


H 


N 


SM 


6 

6 

k 


A 

I) 

3 


k3 


VS 

l) 


sec 

sec 

m 

m 

m 

arc  sec 

arc  sec 

User-defined  constant 
Default  value:  k.(  =  1.74 

User-defined  constant 
Default  value:  ks  =  2.00 

km/sec 

m  (for  definition,  see  algorithm  segment  "(K  I  \\ 

DEPTH  ALGORITHM" 


Algorithm 

N,-  =  ''7  [«,„  -  3  <N,  +  N„)] 


(Mil  1 


13 


NSWC  TR  82-91 


(502) 
(503 1 

( 504) 

Output  Data 

N  .  Geoid  elevation  “observed"  at  the 

seamount  eenter  to  he  matched  tu  the 

calculated  geoid  elevation  during 

parameter  adjustment  m 

Bs  Estimate  of  seamount  half  width  at 

the  base  m 


K  = 


k5  -  h, 

2000 


il)  -  30001 


Bs  =  500  "  "h  -  U  ,Kt> 


=  (1.51  +  05) 


i6ai  +  ,6bi 


i  SO  x  3 <>00 


Initial  estimate  of  seamount  slope 

angle  deg 


OCEAN  DEPTH  ALGORITHM 

l  or  each  data  wind  >w  successfully  processed,  the  ocean  depth.  D.  at  the  seamount  perimeter 
is  required.  Suitable  ocean  depth  data  (Reference  25  in  NSWC  TR  81-200)1  are  available  as 
averages  over  one-by-one  degree  surface  area  elements,  as  a  permanent  data  file  on  the  local 
computer  system.  If  necessary,  convert  D  to  meters. 


ESTIMATION  OE  PEAK  DEPTH 
Subroutine  “DN”  -  Input  Data 


D 

Ocean  depth  at  sea  mount  perimeter 

m 

ds  or  d 

Depth  oi  seamount  peak 

m 

Seamount  slope  angle 

ih'i! 

"s 

Seamount  height 

111 

1 

Crustal  thickness 

m 
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BR  Halt  width  ot'  seamount  root  as  base  m 

ltK  Depth,  below  crust,  of  seamount  root  in 

p^  Seamount  density  gr  in  ' 

p^  Water  density  gr  nr 

pR  Root  density  gi  in 1 

PM  Mantle  density  gr  nr' 

(i  g  Ratio  ot  Newton’s  constant  to  surface  gravity  nr  gi¬ 


lt'  not  otherwise  specified,  assume 

Ps  =  2. bO  l  +t)(i  gr  in3 
p^  =  1.03  l.+0(i  gr/in3 

pR  =  2.05  f+06  gr/ni3 

pM  =  3.40  L;+06  gr,  m3 

Ci,/g  =  0.08024  1.-14  nr.gr 

Subroutine  "1)N”  -  Algorithm 

Reference  is  made  to  l  igure  2.  Perform  the  following  calculations. 


Qs  =  tim 

("0!  i 

=  d,Hs 

i  "02) 

ftR  “  ^R'^R 

I  703  l 

^  =  (1)  +  f  +  1IR  )  Hr 

t  “04) 

evaluate  1  ;U(«S.  ds  1 

(Subroutine  “I  U") 

t  “05 ) 

evaluate  11  (qr  .  0R  ) 

(Subroutine  “IT") 

i  “Oto 
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AXIS  OF  ROTATIONAL 
SYMMETRY 
I 


FIGURE  2.  GEOMETRY  OF  SEAMOUNT  MODEL 


1(. 
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estimation  of  Peak  Depth  for  Isostatieallv  Compensated  Seamount  Input  Data 


t ,  g 

nr  gr 

Ps 

gr  Ill3 

fJ\\ 

gr  in3 

^  K 

gr  m3 

gr  m3 

1) 

m 

T 

m 

vsu 

deg 

«s 

in 

\ 

in 

If  not  otherwise  specified,  assume  for  (i  g  and  the  four  densities  the  nominal  values  listed  under 
"Subroutine  'l)V  Input  Data." 


estimation  of  Peak  Depth  for  Isostatieallv  Compensated  Seamount  Removal  of 
Ill-Conditioned  Cases 


Reference  is  made  to  I  igurc  2.  execute  the  following  test.  This  test  assumes  that  the  sea¬ 
mount  is  as  large  as  the  estimated  seamount  half  width  at  the  base  and  the  local  ocean  depth 

permit  I  seamount  peak  to  be  near  surface  . d*  -  It)  nit.  thus  producing  the  largest 

maximal  gcoid  elevation  D\*  consistent  with  the  particular  Bs  and  I)  values.  If  1)N*  is  smaller 
than  the  observed  maximal  gcoid  elevation  N  . .  the  next  program  segment  will  be  unable  tv' 
perlorm  the  depth  estimation  unless  certain  corrective  measures  are  taken  as  prescribed  below. 


=  ISs 


vl*  =  It)  111 


Ts:  ^  l;in 


-  I 


D-d* 


("i  :• 


t  i t 


(  "141 


IS 


\S\U  IUS201 


I >  i  0  n: 


It* 

it* 

s 

1 

117  - 

Ps  ‘  Pv 

-  II* 

<  " 

K 

D\’ 

-  1  >\t  1).  d'.  v-x.  H*.  1 

i  \>,r .  li*  • 

f  “ 

i  >  \  - 

s  r  s  s 

•  n  : 

k  k 

it  \  os.  proceed  uim  the  next  computer  program  seement.  Otherwise.  prim  out  "(  Mill'S  li  I  - 
U)M)U  IOM  D  (.  \S1  perform  the  assignment  N(  -  DN*.  uiki  consider  the  iust-calce  ,itut  puranu  tern 
marked  witli  a  to  he  the  estimated  seamount  parameters. 


Estimation  of  Peak  Depth  tor  Isostatiealh  compensated 
Seamount  Ueorilhm 

Pe.,'otm  the  following  ealetilalions. 


d  =  D  -  If  tan  y\. 


St  i 


'  \  es.  continue  1 1  no.  eo  to  (  7  23 1. 


d  =10  m 


1)  -  u 


(  ~  24  i 


»R  ^  “s 


"m- - Hsu 

'\l  '  'JK 

l)\|  =  DN").u|.^(,.lls,,.  I.UK.  1 1 K  , , 1 


<  :  ) 


^S  I  "  0  s  X  On 


<  2  M 


I'l 


NS\NC  TK  K2-9! 


‘i  "  l)  »s  tjn  ys  i 


"s,  -l>~ 


i  _'•>  i 

i  ~30i 


«r  ^  «S 


I  -'ll 


PS  ■  *\v 

*R I  "si 

^  R 


(  .'.I 


*>N.  =l>NH).d1.^rllsr  I- Hr- Hk 


("33  l 


Reference  is  made  to  figure  3.  Using  the  above  parameters  with  indices  "0”  and  "I"  as 
initial  values,  enter  now  the  following  iterative  loop: 


vso :  ~  y 


\  " 

■ — ■ —  -  | 

Si  i  vv  _  I  iv  yS  i  *  ! 


U'Vi  -  1A 


("34) 


d„>  =  l>  -  Bs  tan  v-Si+: 
di+:  -  U  ? 

If  ves.  continue.  If  no.  no  to  ("3M. 


t  "35 1 


(  ao  i 


d  ,  =  1 0  111 
!  ♦  2 


(  > 


=  tan 


- 1 


1)  -  d 


( "3M 


^Si+  2  l5s  ' 1,11  S"s i+ : 


bR  Bs 


t"3lM 
( "40 1 


"  *\v 

Hr,  +  ;  =  - - -  HSi+1 

Pm  - 


(741 ) 


l)n(l).  d|+,.  ^Sl+i-  ^r  •  **r  I* :  * 


(742) 
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O.UUOUl  111 


:  ,  -  V  -  D\  , 

1  ♦  -  i  l  * 


Hi, URL  3.  ILLUSTRATION  OF  1TLRATION 

ll  yes.  slop  the  iteration.  Call  the  parameters  associated  with  the  index  i  +  2  the  estimated 
seamount  parameters  in  particular. 

JS  = 

=  ^si+: 

It  no.  reenter  the  iterative  loop  with  the  “i  +  I”  and  "i  +  2”  parameters  as  the  initial  values. 


estimation  of  Peak  Depth  for  lsostatically  Compensated 
Seamount  Output  Data 


dS  =  l‘s,+  s 


estimated  peak  depth  in  m 


1  Estimated  slope  anjile  in  dee 


\S\U  TR  <S2-yi 


"s  =  "s,.: 


II 


k  i '  ' 


I  stimateu  seamount  width  at  h.ise  m  m 
I  stnn.iled  seamount  lieiuht  m  m 
I  stimateu  width  ot  root  in  m 
1  'tunated  lieiuht  of  root  in  m 


in  ease  the  paiameter  estimation  is  terminated  at  I  quation  assian  values  to  the  estimated 

seamount  parameters  as  specified  in  that  latter  ei|uation. 


Lstimation  of  Peak  Depth  for  the  Generally  Compensated 
Seamount  Input  Data 


G  a 

nr  irr 

ur  in3 

ur  m3 

ur  m3 

ur  m3 

D 

m 

I 

tn 

^so 

deg 

Ill 

111 

sk 

"k 

Ill 

*)  "> 
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It  not  otherwise  specified,  assume  for  G  g  and  the  four  densities  the  nominal  values  listed  under 
“Subroutine  'DY  Input  Data." 


Lstimation  of  Peak  Depth  for  Generally  Compensated 
Seamount  Removal  of  Ill-Conditioned  Cases 

Reference  is  made  to  f  igure  2.  f.xecute  the  following  test.  This  test  assumes  that  the  sea¬ 
mount  is  as  large  as  the  estimated  seamount  half  width  at  the  base  and  the  local  ocean  depth 

permit  (seamount  peak  to  be  near  surface  . d*  =  10  mi.  thus  producing  the  hugest 

maximal  geoid  elevation  1)N*  consistent  with  the  particular  Bs  and  D  values.  If  DN*  is  smaller 
than  the  observed  maximal  geoid  elevation  N(. .  the  next  program  segment  will  be  unable  to 
perform  the  depth  estimation  unless  certain  corrective  measures  are  taken  as  prescribed  below. 


II 

»s 

(74b) 

Js  = 

10  m 

i  '47) 

- 

tan-1  ^ 

D  -  d*\ 

(■'481 

^s  ~ 

Bs*  / 

If*  = 

I)  -  10 

m 

( 749) 

B*  =  sk  x  Bs  (750) 

=  »r  <751' 

DN*  =  l)N(l).  d*.  y-*.  II*.  I  .  B*  .  H*)  (752) 

DN*  ^  N(. .'  ("55 1 


If  yes.  proceed  with  the  next  computer  program  segment.  Otherwise,  print  out  “CAUTION: 
ILL-CONDITIONf D  CASf."  perform  the  assignment  N(.  =  DN*.  and  consider  the  just  calculated 
parameters  marked  with  a  “*"  to  be  the  estimated  seamount  parameters. 
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Estimation  of  Peak  Depth  for  the  Generally  Compensated 
Seamount  -  Algorithm 


Perform  the  following  calculations. 


d(1  =  D-  Bs  tan^S() 

(754) 

d»  <  0  ? 

(755) 

If  yes.  continue.  If  no,  go  to  (757). 


d(|  =  1  0  Ill 

*so  =  tan" 


(756) 


(757) 


Hso  =  0“  Jo 

(758) 

BR=skxBs 

(759) 

I)N0  =  DN(D.d0^s0.HS0,T.  Br.Hk) 

(760) 

*si  =  0-8x*so 

(761) 

d,  =  0-  Bs  tan  i?s  j 

(762) 

IIS1  =l)-d, 

(765) 

BR  =  sk  x  Bs 

(764) 

DN,  =  DN(D,d,,^sl,  HS1.T,  BR.  HR) 

( 765 ) 

Using  the  above  parameters  with  indices  “0”  and  “1"  as  initial  values,  enter  now  the  following 
iterative  loop: 


Nc  -  DSL 

^Si+2  ^Si  +  {yfvj  ^  _  [jfvj  ^Si+1  ^Si' 


(766) 
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dt,  =  1)  -  Bs  tan  <'<’') 

d.„,  ^  0  ?  <76*l 

i+  - 


If  yes.  continue.  If  no.  go  to  (770). 


d  ,  =  10  111 
i+  z 


<7<)9) 


^si+:  =  tun 


(770) 


"si*:  lan  ^si+2  (li 

Br  =  sk  x  Bs  (772i 

DN1+:  =  »N(D.  d,+  r  *Slt,.  llSj+2.  T.  Br  .  Hr)  (773. 

e  =  0.00001  Ill  ('74) 

ej+2  =  ,NC  -  i)Nj+,|  ('75) 

ej+,  <  e  ?  ('761 


If  yes.  stop  the  iteration.  Call  the  parameters  associated  with  the  index  i  +  2  the  estimated 


seamount  parameters.  In  particular. 


lls 

^Si+2 

^s 

^Sl+2 

If  no.  reenter  the  iterative  loop  with  the  “i  + 


1"  and  “i  +  2”  parameters  as  the  initial  values. 


Estimation  of  Peak  Depth  for  the  Generally  Compensated 
Seamount  Output  Data 


si* : 


Estimated  peak  depth  in  m 
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v's  ysi+: 
»s 


»R  =  '  HS 

"k 


l  st ii lulled  slope  angle  in  deg 
I  stimated  seamotinl  width  at  base  in  m 
I  stimated  seamount  height  in  in 
1  stimated  width  of  root  in  in 
Specified  height  of  root  in  m 


In  ease  the  parameter  estimation  is  terminated  at  hquation  753.  assign  values  to  the  estimated 
seamount  parameters  as  specified  in  that  latter  equation. 


ANALYSIS  OF  DATA  D ISPLRSION 


Input  Data 


Obtain  the  following  from  the  SIASAT-A  data  track  and  or  data  track  header  (Tables  5  and 
4  of  I  K-8 1-200). 1 

«N  Calculated  average  of  all  oN  (Word  (n  fable  2 1 

on  data  track  under  consideration.  This  is  the 

average  geoid  height  confidence  bound.  m 

Calculated  average  of  all  o§  (Word  7.  1  able  2) 

on  data  track  under  consideration.  This  is  the 

average  vertical  deflection  confidence  bound.  arc  sec 

vs  Satellite  subtrack  velocity  m  sec 

At  Time  interval  between  data  points  see 

further,  have  available 

Larth's  "radius"  ( K r  =  6378000  m>  m 

Ci/g.  ps.  .  pK .  pM  .  as  specified  above. 

1)  Ocean  depth  m 
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Al) 

r 

Al 

o 

N(. 

Bs 

sk 

»K 


Uncertainly  in  ocean  depth  (user  input I. 

Normally,  use  the  value  of  I  O'  of  1).  m 

C  rustal  thickness  (user  input)  in 

Uncertainty  in  crustal  thickness  (user  input)  in 

Initial  estimate  of  slope  angle  deg 

“Observed”  geoid  elevation  above  seamount  center  m 

S cuniouiil  half  width  at  base  m 


(ieometry  factor  for  generally  compensated  cases 
"Height”  of  root  for  generally  compensated  cases  m 


Error  Estimates 

^'s\i  %  ^ 

(SOI  ) 

(S02) 

A5  *  a5 

(Mill 

At A  %  At(j  =  At 

(S05  I 

A^..  *  (1.51  +  05)  AS  deg 

s  ISO  \  3oOO 

1  N0o  1 

ANa  ^  ^  y1\ 

(St)") 

ANC.  *  -  0.5  N;. 

(SOM 

vs  At  ISO 

AXsm  "  ^sm  "  K(  „  lkt 

i  SO1') 

A  IP  %  14  vs  At  1000  m 

(S10) 
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Data  Dispersion  tor  Isostatieally  Compensated  Seamounts 


^isostat.  cornp.^D  +  At) 

i  s  o  $  t  a  t . 

c  o  m  p  . 

(8121 

^  isos  tat.  I’omp.V  +  ^T 

^  i  s  o  s  t  a  t . 

c  o  m  p . 

(813) 

(AJWS  * 

^isostat.  co  m  p .  ^  ^  +  A  \p  ^ 

-  d 

iso  st  a  t . 

c  o  m  p  . 

(814) 

,A^»S  * 

^  isostat.  comp.^B^  +  AB^ 

-  d 

i  s  o  s  t  a  t . 

C  O  111  P  . 

(815) 

^jsostat.  c o  m  p .  ^  +  A  N 

-  d.  ,  , 

i  s  o  s  t  a  t . 

c  o  m  p . 

( 8  1 0> ) 

:«s 
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Data  Dispersion  Analysis  for  Generally  Compensated  Seamounts 

Symboh/e  the  computer  routine  that  performs  the  peak  depth  and  slope  estimation  lor 
generally  compensated  seamounts  by 


,  =  v(D. 

.!  e  d  oral  co  nt  [' . 

aleulate 

r.  ^s.  bs.  \  » 

<  s  r  i 

t  Adi  „ 

^  (d  .  .. 

'  d 

1  N  1  S  1 

a  D 

general  comp.  1>  +  ^  1) 

general  comp. 

<  A  J  >  . 

-  Ui  .  )r 

-  d 

«>!-’> 

A  1 

general,  comp.  1  >  A  1 

;  e  iu‘  r  j  I  v’  *  *  m  {’ . 

t  Ad) 

^  (d  .  1  . 

-  d 

i.s2ui 

•iv's 

o  cm>  c  r  j  1  comp.  v ^  y- 

_•  o  n  o  r  a  l  *.  i «  m  p  . 

t  Ad  I  „ 

^  ( d  ,  l , ,  , , 

-  d 

i  s ;  1 1 

A  1(  s 

o.  c  n  c  r  a  1  c  o  m  p  .  B  M  B  s 

:  n  f  i  4  1  s-'  m  p 

(rid),  , 

^  ( d  .  h  . 

-  d 

t  > i 

A  N  . 

uo  ncra  I  co  in  p  .  \  ^  ^  \  ( 

l  0  71  CM  a  1  c  *  Ml)  p 

Additional  Details  Concerning  the  Data  Dispersion  Analysis 

Note  that  the  two  algoritnms  indicated  by  1  qiiations  Ml  and  M"  depend  on  additional, 
unlisted,  parameters.  Not  specifically  mentioned,  for  example,  are  the  densities.  Also,  sh  and  11 
are  omitted  from  Liquation  817.  All  these  quantities  were  disregarded  because  the  necessary 
quantitative  data  concerning  the  structure  of  the  seamount  root  system  and  the  seamount  compo¬ 
sition  were  not  available  with  the  degree  of  accuracy  that  would  make  inclusion  into  the  data 
dispersion  study  meaningful. 

AUTOMATIC  ANALYSIS  01  SATLLLITL  ALTIMLTRY 

To  perform  an  automatic  seamount  survey  on  Sl.ASAT-A  altimetry,  first  specify  the  desired 
data  track  per  device  set  number  and  revolution  number.  Be  aware  that  a  SI  ASA T-A  altimetry 
track  Is  likely  to  consist  of  several  distinct  track  segments,  l  or  each  of  these  track  segments,  the 
present  seamount  detector  will  automatically  perform  the  Matched  1- liter,  the  RoughnC'  Detector, 
the  Seamount  Locator,  and  the  Calculation  of  Initial  Seamount  Parameters.  Lor  each  trad-  segment, 
indicate  the  cases  rejected  and  those  retained  for  further  processing.  1  or  the  cases  ejected  as 
well  as  for  those  retained,  print  out  the  longitude  ASM  and  the  latitude  y?s M  resulting  liom 
liquations  415  and  416. 

Tor  the  valid  detections,  execute  the  estimation  of  Peak  Depth  for  Isostatically  Compensated 
Seamounts  (this  the  associated  Analysis  of  Data  Dispersion.  Then  execute  the  uncompensated 
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(sk  -  2  aiul  llR  -  0.000001  m)  of  the  Estimation  of  Peak  Depth  for  General!;,  Compen¬ 
sated  Seamounts  plus  the  associated  Analysis  of  Data  Dispersion. 

Make  provision  for  the  user  to  execute,  automatically  for  the  entire  data  track  or.  optionally, 
for  selected  individual  detections,  the  Estimation  of  lVak  Depth  for  tienerally  Compensated 
Seamounts  and  the  associated  Analysis  of  Data  Dispersion  for  arbitrarily  specified  values  of  sk 
and  II  . 

Arrange  the  operational  computer  printout  of  the  input  data  and  results  as  indicated  ny  the 
experimental  data  printouts  reproduced  in  Appendixes  A  through  C.  for  which  credit  is  due  to 
John  l  His  of  the  Physical  Sciences  Software  Branch. 
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APPENDIX  A 

I  RI  AL  DATA  LOR  ISOS  I  V  LK  COMPENSATION 


I 

! 

i 


\S\U  I  K  S2-9I 


SEAMOUNT  DETECTOR 
VERSION  EE 


USER  INPUT 


TYPE  OF  COMPENSATION  =  ISOSTATIC 


CRUST 

T  H  IC  KN  ES  S  - 

5000.00  M 

L 

=  70.00  KM 

CELN 

=  .  020  0  M 

0/G 

=  .63G  2400  0E 

-14  H  **  2/ GR 

RHO.M 

=  3400000.0 

GR/ R**3 

RHO  .R 

=  29  50  0C0.  0 

GR/*'**3 

RHO.  S 

=  26QG0GQ.0 

GR/M**3 

RH  0  .W 

=  1030000.0 

GR/w**3 

K  3  = 

1.740 

K5  = 

2.000 

COG  ROUGHNESS  OETECTOR  will  8E  USEO  WITH 
F  (DATA  W I N 00 W  WIOTh  FACTOR)  =  1.30 

REJECT  DATA  SEGMENTS  NOT  INSIDE  THE  AREA  BOUNDED  BY! 
LATITUDE  =  3  C  60 

LGNGITUOE  =  285  330 

**  NOTE  -  ANY  DATA  SEGMENT  WITH  FEWER  THAN 
8u  DATA  POINTS  WILL  BE  REJECTED. 


>>>>>>>>>>  PROCESS  NEW  DATA  SEGMENT 


header  record 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

CELTA  TIME  =  .490024  SEC 

SAT  lLLITE  VEL. =  6.7840  KM/SEC 

NO.  OF  CATA  POINTS  =  141 

*****  DATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKI»  IT. 


vsn<  i  u  i 


>>>>>>>>>>  PROCESS  NtV  OAT  ft  SEGMENT  <<<<<<<<<< 


READER  9  EC  OR  2 


*EV  NUMBER  =  13  7  5 

PhD  =  .2120 


N  =3 

;LL' A  T I-E  =  .490024  SEC 

SATr'LLIT*-  v/ 1 1_ .  =  6.  7660  KM/SEC 

NJ.  OF  C^TA  POINTS  =  156 


*****  DATA  SEGMENT  NCT  IN  DEFINED  AREA,  SKIP  IT. 


>>>>>>>>>> 


PRGCESS  NEW  DATA  SEGMENT  << 


HEADER  record 


RtV  N CM  PER 
Rh  J 
N 


=  1375 
=  .2120 
=  3 


OciLTA  TIME  =  .490024  SEC 

SAT  tLLl T t  VE L.=  6.  7660  KM/SEC 

NJ.  OP  GATA  POINTS  =  6 


*****  not  Enough  data  points  in  segment,  skip  it. 


PROCESS  NEW  DATA  SEGMENT 


header  record 


REV  MUMPER 

PHO 

N 

DELTA  TIME 


=  1375 
=  .2120 
-  3 

=  .490024  SEC 


SATELLITE  VEL .  =  6.  779  0  KM/SEC 

NO.  OF  CAT  A  POINTS  =  1645 


number  of  data  points  in  each  window  =  22 


NUM3ER  OF  DATA  POINTS  IN  EACH  "COG"  DERIVEO  WINOOW  =  28 


NSWCTR  82-91 


WIN  JON 

NO.  1 

FIRST! 

LATITUDE  =  49 

20 

2 

CEG.HIN.SEC 

LONGITUOE=  30b 

54 

36 

CEG  , MIN, SEC 

LAST  1 

LATITUDE  =  49 

38 

18 

CEG,  MIN, SEC 

LONGITUDE:  306 

17 

32 

CEG, MIN, SEC 

***  PFJECTEO,  MIN.  VERT 

.  OEFL. 

IS  LAST  PT 

OF  WINDOW  **♦ 


PROCESS  NEW  DATA  SEGMENT 


HEADER  RECORD 


REV  NUMPcR 
RHO 

N 


1375 

.2120 

3 


DELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7943  KM/SEC 
NO.  OF  OAT  A  POINTS  =  994 


NUM9ER  OF  DATA  POINTS  IN  EACH  WINDOW  =  22 


NUMdER  OF  DATA  POINTS  IN  EACH  “COG"  DERIVED  WINDOW 
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WINDOW 

NO.  1 

FIRST  t 

L  A  T I TU  OE  -  46 

49 

33 

CEG , MI N , SEC 

LO  NGI TUDE  =  304 

45 

49 

OEG , HI N, SEC 

LAST  : 

LATITUDE  =  46 

7 

21 

ceg,hin,sec 

LONGITUDE =  304 

11 

49 

OEG, MIN, SEC 

***  R^JECTEu,  MIN  •  VERT 

.  OEFL. 

IS  FIRST  PT.  OF  WINDOW  *** 

WINDOW 

NO.  2 

FIRST* 

LATITUQE  =  39 

16 

57 

CEG , MI N,SEC 

L  0  NG  I  7  UCE  =  2  99 

21 

24 

geg,min,sec 

LAST  : 

LATITUDE  =  33 

33 

28 

l£G,MIN,SEC 

LONG  I TUD£=  293 

54 

0 

CEG, MIN, SEC 

ACCEPTE 

D  for  further  analysis 

W IN  DOW 

NG.  3 

FI»ST  : 

LATITUOE  =  32 

45 

59 

OEG, MIN, SEC 

LONG  I  T UCE  =  295 

32 

21 

CEG, MIN, SEC 

LAST  : 

LATITUDE  =  32 

1 

50 

CEG, MIN, SEC 

L  G  NG I T  UOE=  295 

s 

35 
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RESULTS  OF  ANALYSIS  ON  HINOCM  NO.  2 


GEOPHYSICAL  CONSTANTS* 


G/G 

- 

.68  0  240 0 0E- 14 

M**2/GR 

RHO.S 

- 

2600000.0 

G«/H**3 

RHQ.W 

- 

1 0 300  0 u  ••  0 

GR/H**  3 

RHO.R 

- 

2950000.0 

GR /M**  3 

RHO.N 

- 

3400000.0 

GR/M**3 

OCEAN 

OEPTH  ANO  CRUST 

THICKNESS  1 

C 

=  4900.0000 

M 

T 

=  5000.0000 

H 

RESULTS 

FR 

OM  SEAMOUNT 

LOCATOR* 

T  .  A 

145. 74795 

SEC 

T.9 

r 

151 . 65393 

SEC 

T  .SM 

= 

148.80395 

SEC 

N.  A 

= 

-35.552558 

N 

N.  9 

= 

-35.932628 

M 

N.SM 

- 

-34.  796667 

M 

OELTA.A 

= 

-11. 436224 

ARC  SEC 

OELTA  .9 

= 

15.777345 

ARC  SEC 

LONG.SM 

= 

299.14165 

OEG 

299  8  30 

OEG, MIN, SEC 

L AT . SM 

= 

38.942604 

OEG 

38  56  33 

OEG, MIN, SEC 

INITIAL  ESTIMATES* 

fi.S  = 

41153.715 

M 

N.C  = 

1.  49  7  744  8 

M 

PHI. SO  = 

9. 8951 328 

OEG 

ESTIMATE 0 

SEAMOUNT  PARAMETERS* 

ISOSTAT IC 

COMPENSATION 

O.S  = 

109  3.2903 

M 

PHI.S  = 

5.2848116 

OEG 

8.S  = 

41153.  715 

M 

H  .  S  = 

3806.7117 

M 

8.R  = 

41153. 715 

M 

H.R  = 

13291.194 

M 
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Odra  DISPERSION  analysis  on  wincow  NO.  2 


A  VG  . 

C £ 0 1 0  HT. 

= 

«  4124 

37 

31E-0 

1  H 

A  VG  . 

VEkT 

.  OEF 

L  •  - 

.  3313 

54 

06 

ARC  SEC 

SATELLITE 

VEL. 

=  6 

.  78400 

00 

KM/SEC 

TIME 

INTERVAL 

= 

• 

49 

00  24 

SEC 

EARTH 

RAO 

I  US  = 

637  3  0 

00 

.0  H 

UNPER 

TURB 

EO  IN 

PUT  PARAMET 

ER 

St 

C 

4900. 

0000 

4 

T 

- 

5000. 

0000 

M 

PHI  .S 

0  = 

9.995 

1328 

DEG 

8 .  S 

- 

41153 

.  715 

M 

N.C 

- 

1.497 

7448 

M 

ERROR 

EST 

IMATESI 

DELTA 

0 

490. 

00000 

H 

OELTA 

r 

- 

500  . 

00000 

M 

GELTA 

PHI 

.  so  = 

.60457773 

OEG 

OELTA 

3.S 

= 

4654 

.  0519 

H 

OELTA 

N.C 

= 

-.449 

32343 

M 

RESULTS 


parameter  new 

CHANGED  VALUE 

D.  S 

(Mi 

PHI.S 

(OEG) 

H.S  8.R 
(M)  (N) 

H.R 

(M) 

D 

5390 .000 

1491. 

5.412 

3899.  41154. 

13604 

T 

5500.000 

1171. 

5.177 

3729.  41154. 

13009 

PHI. SO 

10. 500 

1093  . 

5  .285 

3807.  41154. 

13281 

B.S 

45807.767 

1239. 

4.569 

3661.  45808. 

12772 

N.C 

1.0  48 

1919. 

4.143 

2981.  41154. 

10400 

CHANGE 
IN  O.S 

(H) 

397.50 

78.154 

0. 

146.09 

825.81 
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RESULTS  OF  ANALYSIS  ON  WINOCW  NC.  3 


*****  WARNING  -  ILL-CONCITICNEO  CASE 


GEOPHYSICAL  constants* 


G/G 

- 

.  68  0  240  00E- 14 

M**2/GR 

RHO.S 

= 

2600000.0 

GR/M**  3 

RHO.W 

1030000.0 

GR/M**3 

RHO.R 

= 

2950000.0 

GR/M**  3 

RHO.  M 

3400000.0 

GR/M  **  3 

OCEAN 

depth  and  crust 

THICKNESS* 

O 

=  4000.0000 

M 

T 

=  5000.0000 

M 

RESULTS 

FROM  SEAMOUNT 

LOCATOR! 

T.A 

z. 

265.53176 

SEC 

T  .  8 

r 

270 . 70577 

SEC 

T  .SM 

r 

268.05602 

SEC 

N.  A 

= 

>42.206865 

H 

N  .  8 

- 

-42. 330503 

H 

N.SM 

-40. 359934 

M 

OtLTA.A 

-33.073740 

ARC  SEC 

DELTA. 8 

r 

33.362061 

ARC  SEC 

LONG. SM 

- 

295.31511 

OEG 

295  19  5 

OEG, MIN, SEC 

L AT . SM 

= 

32.357063 

OEG 

32  21  25 

OEG, MIN, SEC 

INITIAL  ESTIMATES! 

B.S  =  32818.998 

M 

N.C  =  1. 52  329  42 

M 

PHI. SO  =  24.156746 

OEG 

ESTIMATED  SEAMOUNT  PARAMETERS* 
ISOSTATIC  COMPENSATION 

O.S  =  10.000000  M 

PHI.S  =  6.951 7687  OEG 

8.S  =  328 18 . 995  M 

H.S  =  3990.0000  N 

B.R  =  328 18.998  M 

H.R  =  159  20.66  7  N 


59 


NSWCTR  82-91 


OA  T  A  01 SPE RS ION  ANALYSIS  ON  WINOOW  NO.  3 


AVG.  GEOIO  HT.  =  .41243731E- 01  M 
AVG.  VERT.  OEFL.  =  .33135406  ARC  SEC 


SATELLITE  VEL.  =  6.7840000  KM/SEC 

TIME  INTERVAL  =  .490024  SEC 


EARTH  RAOIUS 


6379000  .0  M 


UNPERTURBED  INPUT  PARAMETERS! 


D  = 

4000.  0000 

M 

T  = 

5000. 0000 

M 

PHI . S  0= 

24.  156746 

OEG 

8.S  = 

32818.998 

N 

N.C  = 

3.2446749 

M 

ERROR 

ESTIMATES! 

DELTA 

D 

400.00000 

M 

DELTA 

T  = 

500.00000 

M 

OELTA 

PHr.so= 

.60  45  7773 

OEG 

DELTA 

0.S  = 

4654.0519 

H 

DELTA 

N  .C  = 

-.97346246 

M 

RESULTS 


PARAMETER  NEW  O.S  PHI.S 

CHANGEO  VALUE  (M)  IOEG) 

H.S  B.R 
<M>  IN) 

H.R 

CM) 

CHANGE 
IN  O.S 
IM) 

0 

4400.000  10.  7.619 

WAPNING  -  ILL- CONDITIONED 

4390.  32819. 
CASE  ***** 

15316. 

0. 

T 

5500.000  10.  6.932 

WARNING  -  ILL-CONOITI ONED 

3990.  32819. 
CASE  ***** 

13921. 

0. 

PHI • SO 

24.761  10.  6.932 

WARNING  -  ILu-CONOITIONEO 

3990.  32819. 
CASE  ***** 

13921. 

0. 

6 .  S 

37473.050  10.  6.078 

WARNING  -  ILL-CONOITI CNEO 

3990.  37473. 
CASE  ***** 

13921. 

o. 

N.C 

♦  * 

2.271  10.  6.932 

WARNING  -  ILL-CONOITI CNEO 

3990.  32819. 
CASE  ***** 

13921. 

0. 
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»>»»»>>>>>  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


HEADER  RECORD 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  TIME  =  .490024  SEC 

SATELLITE  V£L.=  6.7800  KM/SEC 

NO.  OF  DATA  POINTS  =  226 

*****  data  segment  not  in  oefineo  area,  skip  it. 


>>>>>>>>>>  process  new  data  segment  <<<<<<<<<< 


HEADER  RECORD 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

DELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7480  KM/SEC 

NO.  OF  DATA  POINTS  =  3048 

*****  OATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT. 


>>>>>>>>>>  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


header  recoro 


REV  NUM8ER  =  1375 

RHO  =  .2120 

N  =3 

DELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7800  KM/SEC 

NO.  OF  DATA  POINTS  =  2206 

*****  DATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT. 


N8ZSE50  ////  END  OF  LIST  //// 
NBZSE50  ////  END  OF  LIST  //// 


41 


NSWC  TR  82-91 


APPENDIX  B 

TRIAL  DATA  FOR  A  TYPICAL  CASE 
OF  GENERAL  COMPENSATION 
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SEAMOUNT  OETECTOR 
VERSION  EE 


USER  INPUT 


TYPE 

OF  COMPENSATION 

=  GENERAL 

CRUST 

THICKNESS  = 

5000.00  M 

SK  = 

2. CO 

HR  = 

3200. 0000 

M 

L 

=  70.00  KM 

OELN 

=  .  0  20  0  M 

G/G 

=  .  680  2  400  OE 

-14  M  **2/GR 

RHO  .M 

=  3400000.0 

GR/F**3 

RHO.R 

=  2950000.0 

GR/W**3 

RHO.S 

=  2600000.0 

GR/M**3 

RHO.W 

=  10  30000.0 

GR/F**3 

K3  = 

1.740 

K  5  = 

2.000 

COG  ROUGHNESS  OETECTOR  WILL  BE  USEO  WITH 
F  (DATA  WlNOOW  WIOTH  FACTOR!  =  1.30 

REJECT  OATA  SEGMENTS  NOT  INSIOE  THE  AREA  BOUNDED  BY« 
LATITUDE  =  30  60 

LONGITUOE  =  285  330 

**  NOTE  -  ANY  OATA  SEGMENT  H TH  FEWER  THAN 
80  DATA  POINTS  WILL  BE  REJECTED. 


»>>>>»»>  PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


MEAOER  RECORO 


REV  NUMBER 

1375 

RHO 

.2120 

N 

3 

DELTA  TIME  = 

.490024 

SEC 

S A  TELL  I Tt  VEL . = 

6. 7840 

KM /SEC 

NO.  OF  OATA  POINTS  =  141 

*****  OATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT. 
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I 

>>>>>>>>>>  process  new  oata  segment  <<<<<<<<<< 

j 

HEADER  RECORD 


REV  NUMBER  = 

1375 

RHO  = 

.2120 

N  = 

3 

DELTA  TIME 

. 49002 

4 

SEC 

SATELLITE  VEL . = 

6.766 

0 

KM/SEC 

NO.  OF  CATA  POINTS  =  1 

56 

*****  OATA  SEGMENT  NOT 

IN 

DEFINED 

>>>>>>>>>>  PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 

HEADER  RECORD 


REV  NUMBER  =  1375 

HO  =  .2120 

=  3 

DELTA  TI*1E  =  .  49  00  2  A  SEC 

SATELLITE  VEL.=  6.7660  KM/SEC 

NO.  OF  OAT  A  POINTS  =  6 

*****  NOT  ENOUGH  DATA  POINTS  IN  SEGMENT,  SKIP  IT. 

>>>>>>>>>>  PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 

HEADER  RECORD 


REV  NUMBER 

1375 

RHO 

= 

.2120 

N 

= 

3 

DELTA  TIME 

.490024 

SEC 

SATELLITE 

VEL  •  - 

6. 7790 

KM /SEC 

NO.  OF  DAT 

A  POI 

NTS  =  1645 

NUMBER  OF 

DATA 

POINTS  IN 

EACH  WINDOW  = 

22 

NUMBER  OF 

OATA 

POINTS  IN 

EACH  '*COG"  OER 

IVEO 

=  29 
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WINDOW 

NO .  1 

FIRST  t 

LATITUDE  =  49 

20 

2 

OEG, 

MIN, SEC 

LONGITUDE =  306 

54 

36 

DEG  , 

MIN, SEC 

LAST  t 

LATITUDE  =  48 

38 

18 

OEG, 

MIN, SEC 

LONG  IT  UOE=  30b 

17 

32 

OEG, 

MIN, SEC 

***  REJtCTEC,  MIN.  VERT.  OEFL.  IS  LAST  PT.  OF  WINDOW  *** 


>>>>>>>>>>  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


HEAOES  RECORD 


REV  NUMBER 

RHO 

N 

DELTA  TIME 


=  1375 
=  .2120 
=  3 

=  .490024 


SEC 


SATELLITE  VEl.-  6.7840  KM/SEC 
NO.  OF  DATA  POINTS  =  994 


NUMBER  OF  DATA  POINTS  IN  EACH  WINOOW  =  22 

NUMBEP  OF  DATA  POINTS  IN  EACH  -COG"  OERIVEO  WINDOW  =  26 
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w : n jo* 

NO .  1 

FWSTl 

LATITUDE  =  46 

49 

38 

OEG, MI N,  SEC 

L0NGITU3E=  304 

45 

49 

0EG.MIN,SEC 

LAST  } 

LATITUDE  =  46 

7 

21 

DEG  »  MIN, SEC 

LONG  I  T  U  ?E  =  3  04 

11 

49 

DEG, MIN, SEC 

***  RLJECTEO,  min.  tf  EkT 

.  OEFL. 

IS  FIRST  PT.  OF  M I NOOW  *** 

WINDOW 

NO.  2 

F  I  R  ;  T  l 

LATITUDE  =  39 

ie 

57 

oeg,min,sec 

LONGITUDE=  299 

21 

24 

oeg,min,sec 

Las"  « 

LATITUDE  =  38 

33 

28 

DEG, MIN, SEC 

LONGITUJE=  298 

54 

0 

OEG , MI N, SEC 

ACCEPTS 

C  FOP  FURTHER  ANALYSIS 

WINDOW 

NO.  3 

FIRST! 

LATITUDE  =  32 

45 

59 

DEG, MIn.SEC 

LO  NG I T  UOE=  295 

32 

21 

ceg,min,sec 

L  a  0  r  : 

LATITUDE  =  32 

1 

50 

OEG, MIN, SEC 

LO Nb  IT  UCE=  295 

9 

35 

ceg,min,sec 

ACCtPTEO  FO  FU°THER  ANALYSIS 
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RESULTS  OF  ANALYSIS  ON  WINDOW  NC.  2 


c-lOohvs ical  o  o  n  s  t  a n t s * 

G/G  =  . 08 C24J0GE-14  M**2/G« 

RHO.S  =  2600300.0  GR/M**3 

RHO.W  -  lUSC'JOii.O  GR/M**3 

RHO.K  =  2960000.0  GR/N**3 

R  h  0 • 1  =  34GG0G0.0  3R/M**3 


OCEAN  DEPTH  Arc;  CRUST  THICKNESS* 

0  =  <*900. 0000  M 

T  =  ST  0  0. 0  0  CO  M 


GENERAL  COMPENSATION  OATA* 

SK  =  2.00uQQ0C 

H.R  =  3200. 0000  M 


RESULTS  FROM  SEAMOUNT  LOCATOR* 


T  .  A 

= 

1  45.  74795 

SEC 

T  .8 

= 

151 . 85393 

SEC 

T  .SM 

148. 80395 

SEC 

N.  A 

- 

-35. 552558 

M 

N.8 

r 

-35.  d  02828 

N 

N  .  SM 

= 

-34. 796667 

M 

CELTA. A 

- 

-11. 436224 

ARC  SEC 

DELTA. 9 

r 

15.777345 

ARC  SEC 

LONG. SM 

- 

299.141 65 

OEG 

299  8  30 

OEG, MIN, SEC 

LAT.SM 

- 

38.942604 

DEG 

38  56  33 

OEG, MI Nt  SEC 

initial  fs 

T I  MATES* 

B.S 

41 153. 715 

M 

N  ,C  = 

1.  4977448 

M 

PHI. SO  = 

9. 8951328 

OEG 

ESTIMATEO  seamount  parameters* 
GENERAL  COMPENSATION 

O.S  -  193  2. 49  62  M 

PHI.S  =  4. 1 243482  OEG 

0.S  =  41153.715  M 

H.S  =  2967.5138  M 

B.R  =  8230  7,4*29  M 

H.R  =  3200.0000  M 
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mu  dispersion  analysis  on  wincow  no.  z 


A  VG  .  GECIO  HI.  =  . 41 ?4 37 3 1  E - 0 1  M 

VERT.  OEFL.  =  .<53135406  ARC  SEC 


SATELLITE  vEL-  =  6.7340000  KM/SEC 

TIRE  INTERVAL  =  .  490024  SEC 


EA^TH  RADIUS  = 


6373000.0  M 


u;,b  E  U'ROE  D  INPUT  PARAMETERS* 


G  = 

4900. 0000 

M 

T  - 

5000. 0000 

i 

UN  I  .  SQ  = 

9.3951323 

DEG 

ri  .  S  = 

41153. 715 

M 

N  ,  C 

1.4977443 

M 

ERROR 

ESTIMATES* 

DELTA 

0 

490.00000 

M 

DELTA 

T  = 

500.00000 

M 

DELTA 

PHI .S0= 

.60457773 

DEG 

DELTA 

8  .  S  - 

4654.0519 

M 

•:  E  L  T  A 

N.C  = 

-  .44932343 

M 

RESULTS 


PARAMETER  NEW 

O.S 

PHI.S 

H.S  8.R 

H  .R 

CHANGED 

value 

(M) 

(DEG) 

(M )  IN) 

IM) 

i) 

5390 .000 

2  378. 

4.186 

3012.  82307. 

3200.0 

T 

5500. 000 

1948. 

4.103 

2952.  82307. 

3200.0 

PHI . SC 

10.500 

1932. 

4.124 

2968.  82307. 

3200.0 

B.  S 

45807.767 

2063. 

3.544 

2837.  91616. 

3200.0 

N.C 

1.048 

2296. 

3.620 

2604.  82307. 

3200.0 

CHANGE 
IN  O.S 

(M) 

445.70 

15.330 

0. 

130.09 

363.78 
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RESULTS  CF  ANALYSIS  ON  wINqOH  no.  3 


WARNING  -  ILL-CONOITICNEO  CASE 


geophysical  constants* 


G/G  = 

. 66024000E-14 

M**2/GR 

RHO.S  = 

2600000.0 

GR/M**3 

RHO.W  = 

1030000.0 

OR/M**  3 

RHO.R  = 

2950000.0 

GR/M**3 

RHO.H  = 

3400000.0 

OR/M** 3 

OCEAN  OEPTh 


ano  crust  thickness* 


0 

=  4000.0000 

H 

T 

=  5000.0000 

M 

GENER 

AL  COMPENSATION 

OAT  A* 

SK  = 

2.0000000 

H.R  = 

3200. 0000 

M 

RESULTS 

FROM  SEAMOUNT 

LOCATOR  * 

T  .  A 

=  265.53176 

SEC 

T.B 

=  270.70577 

SEC 

T.SM 

=  2 68.05602 

SEC 

N.  A 

=  -42.206865 

M 

N.  8 

=  -42.330503 

N 

N.SM 

=  -40.359934 

M 

OELTA. A 

=  -33.0  73740 

ARC  SEC 

DELTA. 9 

=  33.362081 

ARC  SEC 

LONG. SM 

=  295.31811 

DEG 

295  19  5 

DEG, MIN,  SEC 

LAT  .SM 

=  32.  357068 

OEG 

32  21  25 

OEG.MIN, SEC 

INITIAL 

ESTIMATES* 

B.S  =  32016.998  M 

N.C  =  2.3992570  M 

PHI. SO  =  29. 1 5C746  DEG 


ESTIMATES  StAMOuNT  PARAMETERS! 
general  COMPENSATION 


D.S 

= 

10. 000000 

M 

PHI.S 

= 

6.931 7687 

DEG 

B.S 

= 

32618.998 

M 

H.S 

399  0.  0C00 

M 

B.R 

= 

65637. 996 

M 

H.R 

= 

3200. 0000 

M 

*  *  * 
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OATA  DISPERSION  ANALYSIS  ON  WINCOW  NO.  3 


AVG.  GECI3  HT.  =  . 4124 37 31 E- 01  M 
AVG.  VERT.  OEFL.  =  .83135406  ARC  SEC 

SATELLITE  VEL.  =  6.7840000  KM/SEC 

TIME  INTERVAL  =  .490024  SEC 


EARTH  RADIUS 


6373000.0  M 


UNPERTURBED  INPUT  PARAMETERS! 


D  = 

4000.  0000 

M 

T  = 

5000. 0000 

M 

PHI.S0= 

2  4.  156746 

OEG 

B.S  = 

32818.998 

M 

N.C  = 

3.2448749 

H 

error 

ESTIMATES » 

DELTA 

0 

400.00000 

M 

DELTA 

T  = 

500.00000 

M 

DELTA 

PHI.S0= 

.60457773 

CEG 

DELTA 

a.s  = 

4654.0519 

M 

DELTA 

N.C 

-.97346246 

M 

RESULTS 


CHANGE 


PARAMETER  NEW  O.S 

PHI.S 

H.S  8.R 

H.R 

IN  D.S 

CHANGED  VALUE  (Ml 

(OEG) 

(Ml  ( M) 

(Ml 

(Ml 

D 

4400.000  10. 

7.619 

4390.  65638. 

3200.0 

0. 

WARNING  -  ILL-CONUITI ONEO 

CASE  ***** 

T 

5500.000  10. 

6  .932 

3990.  65638. 

3200.0 

0. 

WARNING  -  ILL-CON 0ITICN60 

CASE  ***** 

PHI, SO 

24.761  10. 

6.932 

3990.  65638. 

3200.0 

0. 

WARNING  -  ILL-CONOITICNED 

CASE  ***** 

e.s 

37473.050  10. 

6.078 

3990.  74946. 

3200.0 

0. 

*  ♦  ♦ 

WARNING  -  ILL-CONOITIGNEO 

CASE  ***** 

N.C 

2.271  82. 

6.809 

3918.  65638. 

3200.0 

71.523 
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>>>>>>>>>>  PRJCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


HEAjtR  RECORD 


REV  NUMBER  =  1 3  75 

RHO  =  .2120 

N  =3 

DEL  '  A  TI-E  =  .490024  SEC 

SATELLITE  VEL  .  =  6.  7800  KM/SEC 

NO.  OF  OATA  POINTS  =  226 

*****  OATA  SEGMENT  NOT  IN  DEFINEO  AREA*  SKIP  IT. 


>>>>>>>>>>  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


HEACER  PECGRD 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  TIME  =  .490024  SEC 

SATELLITE  VEL»=  6.7481  KM/SEC 

NO.  OF  OATA  POINTS  =  30  48 

*****  DA” A  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT. 


>>>>>>>>>>  PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


HEADER  RECORD 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  TIME  =  .490024  SEC 

SATELLITE  VEl.=  6.7600  KM/SEC 

NO.  OF  CAT  A  POINTS  =  2206 

*****  OATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT, 


*********** 


N8ZSE5N  ////  END  OF  LIST  //// 
N8ZSE5N  ////  END  OF  LIST  //// 
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APPENDIX  C 

TRIAL  DATA  FOR  ABSENCE  OF  ROOT 
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SEAMOUNT  OET  ECTOR 
VERSION  EE 


USER  INPUT 


TYPE  OF  COMPENSATION  =  GENERAL  -  UNCOMPENSATED 
CRUST  THICKNESS  =  5000.00  M 

SK  =  2. CO 

HR  =  .  10  00  00  0  0E- 05  M 

L  =  70.00  KM 

OELN  =  .0200  M 

G/G  =  .6  80  2  40  0  QE-14  M**2/GR 

RHO.M  =  3400000. 0  GR/N**3 

RHO.R  =  2950000.0  GR/M**3 

RHO.S  =  2600000.0  GR/P**3 

RHO.M  =  1030000.0  GR/W**3 

K3  =  1.740 

K5  =  2.000 

COG  ROUGHNESS  OETECT OR  WILL  BE  USEO  WITH 
F  (DATA  WINOOW  WIDTH  FACTOR)  =  1.30 

REJECT  DATA  SEGMENTS  NOT  INSIOE  THE  AREA  BOUNDED  BY* 
LATITUDE  =  30  60 

LONGITUDE  =  285  330 

**  NOTE  -  ANY  DATA  SEGMENT  WITH  FEWER  THAN 
80  OATA  POINTS  WILL  BE  REJECTED. 


>>>>>>>>>>  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


HEAOER  RECCRO 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

DELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7840  KM/SEC 

NO.  OF  OATA  POINTS  =  141 

*****  OATA  SEGMENT  NCT  IN  OEFINED  AREA,  SKIP  IT. 
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»»>>>>»  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


HEADER  RECORD 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  time  =  .490024  SEC 

SATELLITE  y*EL.=  6.7660  KM/SEC 

NO.  OF  DATA  POINTS  =  156 

*****  DATA  SEGMENT  NOT  IN  QEFINEO  AREA,  SKIP  IT. 


»»>>»>>  PROCESS  NEW  oata  SEGMENT  <<<<<<<<<< 


HEAOER  RECORO 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7660  KM/SEC 

NO.  OF  OATA  POINTS  =  0 

*****  NOT  ENOUGH  OATA  POINTS  IN  SEGMENT,  SKIP  IT. 


>>>>>>>>>>  FRQCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


HEAOER  RECORO 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  TIME  =  .490024  SEC 

SATELLITE  YEL.=  6.7790  KM/SEC 

NO.  OF  OATA  POINTS  =  1645 

NUMdER  OF  OATA  POINTS  IN  EACH  WINDOW  =  22 

NUMBER  OF  OATA  POINTS  IN  EACH  “COG**  DERIVED  WINDOW  =  25 
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WlNOOW 

NO.  1 

FIRST! 

LATITUOE  = 

49 

20 

2 

OEG.MXN.SEC 

LONGI T  UO£  = 

306 

54 

36 

DEG. Ml N.SEC 

LAST  * 

LATITUOE  = 

43 

38 

18 

DEG «  MI N. SEC 

LONGI TUOE= 

3  06 

17 

32 

DEG. MIN, SEC 

***  REJECTEO,  MIN.  VERT.  OEFL .  IS  LAST  PT .  OF  WlNOOW  *** 


>>>>>>>>>>  PROCESS  NEW  DATA  SEGMENT  <<<<<<<<<< 


HEACER  record 


REV  NUMBER  =  1375 

RHO  =  •  2120 

N  =3 

DELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7840  KM/SEC 
NO.  OF  OATA  POINTS  =  994 

NUM8ER  OF  DATA  POINTS  IN  EACH  WlNOOW  =  22 

NUM3ER  OF  OATA  POINTS  IN  EACH  "COG”  DERIVED  WINDOW  =  28 
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W  IN30W 

NO.  1 

FIRST l 

LATITUDE  = 

46 

49 

36 

PEG, MIN, SEC 

LONG  IT  UOE  = 

304 

45 

49 

0£G,MIN,SEC 

LAST  t 

LATITUOE  = 

46 

7 

21 

OEG  , MIN, SEC 

LONG  I T UOE= 

3  04 

11 

49 

DEG, MI N, SEC 

***  REJECTED,  MIN. 

YERT 

.  OEFL. 

IS  FIRST  PT.  OF  WINDOW  *** 

WINDOW 

NO.  2 

FIRST  » 

LATITUOE  = 

39 

16 

57 

DEG,  MIN, SEC 

L  ONG IT  U0£  = 

2  99 

21 

24 

DEG, MIN, SEC 

LAS  r  t 

LATITUOE  = 

33 

33 

28 

OEG,  MIN,  SEC 

LONG  I TUOE= 

2  93 

54 

0 

OEG, MIN, SEC 

ACCEPTEO  FOR  FURTHER  ANALYSIS 

WINDOW 

NO.  3 

FIRST l 

LATITUOE  = 

32 

45 

59 

OEG,  MIN,  SEC 

LONGITUDE= 

2  95 

32 

21 

DEG, MIN, SEC 

LAST  « 

LATITUOE  = 

32 

1 

50 

OEG, MIN, SEC 

L0NGITU0E= 

295 

6 

35 

OEG, MIN, SEC 

ACCEPTEO  FOR  FURTHER  ANALYSIS 
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RESULTS  OF  ANALYSIS  ON  HI  NO  CN  NC«  2 


GEOPHY 

SICAL  CONSTANTS! 

G/G 

=  . 68 024000E-14 

M**2/GR 

RHO.S 

=  2600000.0 

GR/N**3 

RHO.W 

=  1030000.0 

GR/M**3 

RHO.R 

=  2950000.0 

GR/N**3 

RHO.M 

=  3400000.0 

GR/M**3 

OCEAN 

OEPTH  AND  CRUST 

THICKNESS! 

0 

=  4900.0000 

M 

r 

=  6000.0000 

M 

general  COMPENSATION  oatai 

SK  =  2.0000000 

H.R  =  .10000000E-05  M 


RESULTS  FROM  SEAMOUNT  LOCATOR! 


T.A  = 

145. 74795 

SEC 

T .  3 

151.35393 

SEC 

T  ,  S  M  * 

148.80395 

SEC 

N.  A  = 

-35.552558 

N 

N .  3  ~ 

-35.  602823 

M 

N.SM  = 

-34. 796667 

M 

OELTA.A  = 

-11. 436224 

ARC  SEC 

DELTA. 8  = 

15,  777345 

ARC  SEC 

LONG. SN  = 

299.14165 

OEG 

299  8  30 

OEG,HlN,  SEC 

LAT.SM  = 

33.  942  604 

OEG 

33  56  33 

OEG, MIN, SEC 

INITIAL  ESTIMATES! 

0.S  = 

41153. 715 

M 

N.C  = 

1.49  77448 

M 

PHI. SO  = 

9.  3951323 

OEG 

ESTIMATED 

SEAMOUNT  PARA 

METERS« 

general 

COMPENSATION 

-  NO  R COT 

O.S 

PHI.S 

e.s 

H.S 

b.r 

H.R 


3568.2760  M 

1.8256144  OEG 

41153.715  M 

1311.7240  M 

32307.429  M 

. 1COOOOOOE-05  1 
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DATA  DISPERSION  analysis  on  nincow  NO.  2 


A  VG.  GEOIO  HT.  =  . 91 24 37 31E- 01  M 
AVG.  VERT.  OEFL.  =  .83135406  ARC  SEC 

SATELLITE  VEL.  =  6.7840000  KM/SEC 

TIME  INTERVAL  =  .490024  SEC 


EARTH  RAOIUS 


6378000.0  H 


UNPERTURBED  INPUT  PARAMETERS* 


0  =  4900.0000 

T  =  5000.0000 

PHI.SG=  9.8951328 
B.S  =  41153.715 

N.C  =  1.4977448 


M 

M 

□  EG 

M 

M 


ERROR  ESTIMATES* 


DELTA  0  = 
DELTA  T  = 
DELTA  PHI  «  SO  * 
DELTA  B.S  = 
DELTA  N.C  = 


490.00000 

500.00000 

.60457773 

4654.0519 

.44932343 


M 

M 

CEG 

M 

M 


RESULTS 


PARAMETER  NEW 
CHANGEO  VALUE 

D.S 

(Ml 

PHI.S 
(OEG ) 

H.S  B.R 

(Ml  (Ml 

H.R 

CM) 

CHANGE 

IN  D.S 

(Ml 

0 

5390.000 

4052. 

1.862 

1338.  82307. 

. 10000E-05 

464.07 

T 

5500.000 

3588. 

1.826 

1312.  82307. 

.  10000E-05 

.541 33E 

PHI .SO 

10.500 

3538. 

1.826 

1312.  82307. 

.1 OOOOE-05 

0. 

B.S 

45807.767 

3740. 

1 .450 

1160.  91616. 

. 1 0000E-35 

151.81 

N.C 

1.048 

3976. 

1.287 

924.  82307. 

. 1 OQOOE-05 

387.39 
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RESULTS 

OF  ANALYSIS  ON 

R I  NO  CM  no. 

GEOPHYSICAL  CONSTANTS! 

G/G  = 

.  6o  0  240  00  £-  14 

M**2/GR 

RHO.S  = 

2600000. 0 

G  R  /  M  *  *  3 

RHO.W  = 

1  0  30  CO  0.0 

GR/M*» 3 

RHO.p  = 

2950000. 0 

GR/M**3 

RHO.M  = 

340  000  3.  0 

GR/M**3 

OCEAN 

DEPTH  ANO  CRUST 

THICKNESS « 

-  4000.0000 

M 

T 

-  5000.0000 

1 

GENERAL  COMPENSATION  OATAl 

SK  =  2.0000000 

H.R  =  .10  00  300  0£-  05  M 


RESULTS  FROM  SEAMOUNT  LOCATOR! 


T  .A 

- 

265. 53176 

SEC 

T.8 

270. 70577 

SEC 

T  .SM 

- 

260.05602 

SEC 

N  .  A 

r 

-42.206665 

M 

N  .  8 

- 

-42. 330503 

M 

N.SM 

- 

-40. 359934 

H 

DELTA.  A 

r 

-33.073740 

ARC  SEC 

OELTA  .  8 

- 

33. 362061 

ARC  SEC 

LONG. SM 

- 

295.31811 

□EG 

295  19  5 

0£G,  MIN,  SEC 

LAT. SM 

= 

32.357068 

OEG 

32  21  25 

0EG»  MINf  SEC 

INITIAL  ESTIMATES! 

B.S  =  32818.998 

M 

N.C  =  3.  2446749 

M 

PHI. SO  =  24.156746 

OEG 

ESTIMATED  SEAMOuNT  PARAMETERS! 


GENERAL 

COMPENSATION 

-  NO 

ROOT 

O.S 

— 

613.88354 

M 

PHI.S 

a 

5.8906771 

OEG 

B.S 

= 

3  26  18.998 

M 

H.  S 

■= 

3  3o  6. 1165 

M 

8 .  R 

r 

65637.996 

M 

H.R 

= 

.1QOOOOOOE-05 

M 

6: 
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DATA  DISPERSION  ANALYSIS  3N  WlNcOW  NO.  3 


A  VG  .  GEOIO  HT.  =  »  41 2437  31E- 01  M 

A V/G .  VERT.  OEFL •  =  .33135406  ARC  SEC 

Satellite  vel.  =  6.7343000  km/sec 

TIME  INTERVAL  =  .490024  SEC 

EARTH  RADIUS  =  6373000.0  M 


UNPERTURBED  INPUT  PARAMETERS* 

O  =  4000.  0000  M 

T  =  5000.0000  1 

PHI.SC=  24.156746  DEG 

B.S  =  32018.990  i 

N.C  =  3.2448749  M 


ERROR  ESTIMATESi 

DELTA  D  =  400.00000  M 

DELTA  T  =  500.00000  M 

DELTA  PHI . SO  =  .60457773  CEG 

CELTA  B.S  =  4654.0519  M 

DELTA  N.C  =  -.97346246  M 


RESULTS 


PARAMETER  NEW 
CHANGEO  VALUE 

G.S 

(Ml 

PHI.S 

(OEGI 

H.S  B.R 

(Ml  <*> 

H.R 

<M> 

CHANGE 

IN  D.S 

-  —  —  -*■  * 

*  ~  — 

0 

4400. 00C 

943. 

6.013 

3457.  65638. 

•  1 0000E-05 

329. 17 

T 

5500 .000 

614. 

5.891 

3396.  65638. 

. 10000E-35 

. 5820  8£ 

PHI. SO 

24. 761 

614. 

5  .891 

3366.  65638. 

.10000E-05 

t. 

B.S 

37473.050 

1055. 

4.494 

2945.  74946. 

. 1000QE-05 

440.70 

N.C 

2.271 

1571. 

4.232 

2429.  65638. 

. 10000E-35 

957.46 
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PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


HEAOER  RECORO 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  I  IKE  =  . 490024  SEC 

SATELLITE  VEL.=  6.7800  KM/SEC 

NO.  OF  DATA  POINTS  =  226 

*****  oata  segment  not  in  defined  area,  skip  it. 


PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


HEAOER  RECORO 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

DELTA  TIME  =  .490024  SEC 

SATELLITE  VEL.=  6.7480  KM/SEC 

NO.  OF  OATA  POINTS  =  3048 

*****  DATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT. 


PROCESS  NEW  OATA  SEGMENT  <<<<<<<<<< 


HEADER  RECORD 


REV  NUMBER  =  1375 

RHO  =  .2120 

N  =3 

OELTA  TIME  =  .490024  SEC 

SATELLITE  VEL . =  6.7300  KM/SEC 

NO.  OF  DATA  POINTS  =  2206 

*****  OATA  SEGMENT  NOT  IN  DEFINED  AREA,  SKIP  IT, 


*********** 

*•*•**•*«** 


N9Z5E5L  ////  END  OF  LIST  //// 
N9ZSE5L  ////  ENO  OF  LIST  //// 
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